Hybrid quantum-classical modeling of quantum dot devices

The design of electrically driven quantum dot devices for quantum optical applications asks for modeling approaches combining classical device physics with quantum mechanics. We connect the well-established fields of semi-classical semiconductor transport theory and the theory of open quantum systems to meet this requirement. By coupling the van Roosbroeck system with a quantum master equation in Lindblad form, we introduce a new hybrid quantum-classical modeling approach, which provides a comprehensive description of quantum dot devices on multiple scales: It enables the calculation of quantum optical figures of merit and the spatially resolved simulation of the current flow in realistic semiconductor device geometries in a unified way. We construct the interface between both theories in such a way, that the resulting hybrid system obeys the fundamental axioms of (non-)equilibrium thermodynamics. We show that our approach guarantees the conservation of charge, consistency with the thermodynamic equilibrium and the second law of thermodynamics. The feasibility of the approach is demonstrated by numerical simulations of an electrically driven single-photon source based on a single quantum dot in the stationary and transient operation regime

Beyond just ``flattening the curve'': Optimal control of epidemics with purely non-pharmaceutical interventions

When effective medical treatment and vaccination are not available, non-pharmaceutical interventions such as social distancing, home quarantine and far-reaching shutdown of public life are the only available strategies to prevent the spread of epidemics. Based on an extended SEIR (susceptible-exposed-infectious-recovered) model and continuous-time optimal control theory, we compute the optimal non-pharmaceutical intervention strategy for the case that a vaccine is never found and complete containment (eradication of the epidemic) is impossible. In this case, the optimal control must meet competing requirements: First, the minimization of disease-related deaths, and, second, the establishment of a sufficient degree of natural immunity at the end of the measures, in order to exclude a second wave. Moreover, the socio-economic costs of the intervention shall be kept at a minimum. The numerically computed optimal control strategy is a single-intervention scenario that goes beyond heuristically motivated interventions and simple "flattening of the curve". Careful analysis of the computed control strategy reveals, however, that the obtained solution is in fact a tightrope walk close to the stability boundary of the system, where socio-economic costs and the risk of a new outbreak must be constantly balanced against one another. The model system is calibrated to reproduce the initial exponential growth phase of the COVID-19 pandemic in Germany

Non-isothermal Scharfetter--Gummel scheme for electro-thermal transport simulation in degenerate semiconductors

Electro-thermal transport phenomena in semiconductors are described by the non-isothermal drift-diffusion system. The equations take a remarkably simple form when assuming the Kelvin formula for the thermopower. We present a novel, non-isothermal generalization of the Scharfetter? Gummel finite volume discretization for degenerate semiconductors obeying Fermi?Dirac statistics, which preserves numerous structural properties of the continuous model on the discrete level. The approach is demonstrated by 2D simulations of a heterojunction bipolar transistor

Mathematical models: A research data category?

Mathematical modeling and simulation (MMS) has now been established as an essential part of the scientific work in many disciplines and application areas. It is common to categorize the involved numerical data and to some extend the corresponding scientific software as research data. Both have their origin in mathematical models. In this contribution we propose a holistic approach to research data in MMS by including the mathematical models and discuss the initial requirements for a conceptual data model for this field

Beyond just “flattening the curve”: Optimal control of epidemics with purely non-pharmaceutical interventions

When effective medical treatment and vaccination are not available, non-pharmaceutical interventions such as social distancing, home quarantine and far-reaching shutdown of public life are the only available strategies to prevent the spread of epidemics. Based on an extended SEIR (susceptible-exposed-infectious-recovered) model and continuous-time optimal control theory, we compute the optimal non-pharmaceutical intervention strategy for the case that a vaccine is never found and complete containment (eradication of the epidemic) is impossible. In this case, the optimal control must meet competing requirements: First, the minimization of disease-related deaths, and, second, the establishment of a sufficient degree of natural immunity at the end of the measures, in order to exclude a second wave. Moreover, the socio-economic costs of the intervention shall be kept at a minimum. The numerically computed optimal control strategy is a single-intervention scenario that goes beyond heuristically motivated interventions and simple “flattening of the curve”. Careful analysis of the computed control strategy reveals, however, that the obtained solution is in fact a tightrope walk close to the stability boundary of the system, where socio-economic costs and the risk of a new outbreak must be constantly balanced against one another. The model system is calibrated to reproduce the initial exponential growth phase of the COVID-19 pandemic in Germany. © 2020, The Author(s)

Discretization scheme for drift-diffusion equations with a generalized Einstein relation

Inspired by organic semiconductor models based on hopping transport introducing Gauss-Fermi integrals a nonlinear generalization of the classical Scharfetter-Gummel scheme is derived for the distribution function F(η)=1/(exp(-η)+γ). This function provides an approximation of the Fermi-Dirac integrals of different order and restricted argument ranges. The scheme requires the solution of a nonlinear equation per edge and continuity equation to calculate the edge currents. In the current formula the density-dependent diffusion enhancement factor, resulting from the generalized Einstein relation, shows up as a weighting factor. Additionally the current modifies the argument of the Bernoulli functions

Make the most of your visual simulation data

There still exists no common standard on how to publish and cite visualisations of simulation data. We introduce the TIB AV-Portal as a sustainable infrastructure for audio-visual data using a combination of digital object identifiers (DOI) and media fragment identifiers (MFID) to cite these data in accordance with scientific standards. The benefits and opportunities of enhancing publications with visual data are illustrated by showing a use case from opto-electronics

Numerical simulation of carrier transport in semiconductor devices at cryogenic temperatures

At cryogenic temperatures the electron-hole plasma in semiconductor materials becomes strongly degenerate, leading to very sharp internal layers, extreme depletion in intrinsic domains and strong nonlinear diffusion. As a result, the numerical simulation of the drift-diffusion system suffers from serious convergence issues using standard methods. We consider a one-dimensional p-i-n diode to illustrate these problems and present a simple temperature-embedding scheme to enable the numerical simulation at cryogenic temperatures. The method is suitable for forward-biased devices as they appear e.g. in optoelectronic applications

Mathematical models: A research data category?

Mathematical modeling and simulation (MMS) has now been established as an essential part of the scientific work in many disciplines and application areas. It is common to categorize the involved numerical data and to some extend the corresponding scientific software as research data. Both have their origin in mathematical models. In this contribution we propose a holistic approach to research data in MMS by including the mathematical models and discuss the initial requirements for a conceptual data model for this field

Convergence of a finite volume scheme to the eigenvalues of a Schrödinger operator

We consider the approximation of a Schrödinger eigenvalue problem arising from the modeling of semiconductor nanostructures by a finite volume method in a bounded domain $OmegasubsetR^d$. In order to prove its convergence, a framework for finite dimensional approximations to inner products in the Sobolev space $H^1_0(Omega)$ is introduced which allows to apply well known results from spectral approximation theory. This approach is used to obtain convergence results for a classical finite volume scheme for isotropic problems based on two point fluxes, and for a finite volume scheme for anisotropic problems based on the consistent reconstruction of nodal fluxes. In both cases, for two- and three-dimensional domains we are able to prove first order convergence of the eigenvalues if the corresponding eigenfunctions belong to $H^2(Omega)$. The construction of admissible meshes for finite volume schemes using the Delaunay-Voronoï method is discussed. As numerical examples, a number of one-, two- and three-dimensional problems relevant to the modeling of semiconductor nanostructures is presented. In order to obtain analytical eigenvalues for these problems, a matching approach is used. To these eigenvalues, and to recently published highly accurate eigenvalues for the Laplacian in the L-shape domain, the results of the implemented numerical method are compared. In general, for piecewise $H^2$ regular eigenfunctions, second order convergence is observed experimentally